Specific heat of the Kelvin modes in low temperature superfluid turbulence 
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It is pointed out that the specific heat of helical vortex line excitations, in low temperature superfluid 
turbulence experiments carried out in helium II, can be of the same order as the specific heat of the 
phononic quasiparticles. The ratio of Kelvin mode and phonon specific heats scales with Lo T~ 5 ^ 2 , 
where Lo represents the smoothed line length per volume within the vortex tangle, such that the 
contribution of the vortex mode specific heat should be observable for Lo = 10 6 • ■ • 10 s cm -2 , and 
at temperatures which are of order 1 • • • 10 mK. 
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The problem of superfluid turbulence has seen a re- 
newed upsurge of interest, largely due to the fact that 
turbulence in superfluids at temperatures of order one 
Kelvin [0-j3|, has proved to be similar to turbulence in 
normal fluids |Q. Of particular interest should be the 
range of lower temperatures, of order a few mK, where 
the normal fluid density is negligible, and a particularly 
interesting behaviour of superfluid turbulence in a very 
pure form should be observable Since the pioneer- 

ing work of Vinen , the dynamics of superfluid turbu- 
lence, especially its generation and decay, is to a large 
extent still not understood. More recent efforts in this 
direction, investigating, e.g., the energy spectrum and re- 
connection events in the tangle, can be found in ■ An 
interesting thermodynamic approach to directly observ- 
able effects of dissipation in low temperature superfluid 
turbulence has been proposed by Samuels and Barenghi 
pot - The effect to be observed, according to this pro- 
posal, consists in the conversion of (incompressional) ki- 
netic energy of the vortex tangle into compressional en- 
ergy, which then heats the sample. In that work, it has 
been assumed that the specific heat of the turbulent he- 
lium II sample is given solely by phononic excitations. 
Here, it will be shown that, at temperatures of the order 
a few mK, the total specific heat of line excitations in the 
vortex tangle is of the same order as the specific heat of 
phonons, integrated over the volume of the tangle. This 
will be true for large, yet experimentally feasible, line 
densities in the tangle. Though vortex lines fill out a 
very small fraction of the volume of the turbulent helium 
sample, the contribution of their oscillations to the spe- 
cific heat should, as will now be shown, be discernible 
from the phonon specific heat at low temperatures. 

The total specific heat of the helium II sample is as- 
sumed to be composed of the specific heat of the elemen- 
tary nontopological excitations which constitute the heat 
bath, living in the three dimensions of the sample, and 
the specific heat of the excitations propagating along the 
topological line vortices, which live in one spatial dimen- 
sion. Below temperatures of a fraction of one Kelvin, 
only phonons contribute to the specific heat of the non- 
topological excitations, the contribution of the gapful ex- 



citations (rotons) being exponentially suppressed. The 
specific heat of sound excitations in a volume V is 
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where c s is the speed of sound. The specific heat of the 
oscillation modes of an isolated vortex ring of radius r$ 
may generally be written in the form 
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The nondimcnsional function / depends on the ratio of 
the temperature and the 'cyclotron' energy of the vor- 
tex core, ksT/huJc, characterizing the separation be- 
tween thermal and core energy physics |pT| . In helium 
II, huj c ~ 1 • • • 10 K, so that the parameter ksT/huic ~ 
1CP 3 • • • 10~ 4 for temperatures in the mK range. The 
cyclotron frequency of the core is defined by uo c — 
YJ (7r(£ c /2) 2 ) , with £ c the core diameter and T = h/m the 
quantum of velocity circulation. For a line of arbitrary 
global shape, 27rro in (§) is replaced by the smoothed 
arc length of vortex line, i.e., the length of line with the 
Kelvin disturbances subtracted This undisturbed, 
smoothed length of line plays the role of "volume" in the 
specific heat of the Kelvin modes. The smoothed line 
density of vortices is designated Lq = Lo[cm~ 2 ], so that 
we may assign 2nro = LqV, as an (average) filament 
length entering the specific heat (|J). The line density 
Lq, around which there are Kelvin fluctuations, is to be 
distinguished from the full line density L, which derives 
from the total (incompressible) kinetic energy, with the 
contribution of the Kelvin waves included. They differ 
by a factor logarithmically dependent on the product of 
the largest Kelvin wave number excited and the average 
vortex element distance I = Lp ' (the smallest possible 
Kelvin wave number is l^ 1 ) H, so that Lq can be about 
an order of magnitude less than L. 

The function / in (^) depends on the spectrum of waves 
on the vortex line. For wavelengths much larger than 
the core diameter, a single ring vortex has the Kelvin 
spectrum 



1 



LOK 



Tn 2 
47rrg 



In [n*/n] 



(3) 



where n is the number of waves on the circumference 
(the mode number), so that the wave number k = n/r . 
The ultraviolet cutoff mode number is for a vortex ring 
ri* = O(27rr /^ c ); for a filament of more general shape, 
the cutoff wave number k c — 0(^ ( T 1 ), and n* is of order 
the total length of a filament, divided by £ c . 

If vortex filaments approach each other to within a 
distance corresponding to the inverse Kelvin wave vector 
under consideration, so that kl < 1, the Kelvin spectrum 
will be changed. In the expression (^) for the specific 
heat of the oscillations presented below, we will assume 
that vortices remain, at least on average, sufficiently well 
separated to retain the validity of the Kelvin spectrum 
(||) for most of the excitations on the filaments. Devia- 
tions from that spectrum are then too small to be ther- 
modynamically relevant for the specific heat of the os- 
cillations, and in particular its temperature dependence. 
Conversely, if there are any measurable deviations from 
the temperature dependence of the specific heat of Kelvin 
waves, we will have an indicator that the Kelvin spec- 
trum has changed due to the presence of a very dense 
vortex tangle. Strictly spoken, these conclusions can only 
claim validity if the vortex system is in thermodynamic 
equilibrium with itself and the surrounding fluid. This 
is the case for a vortex array, generated by a constant 
rotation rate ft of the superfluid, with resulting densi- 
ties Lq = 2fi/T = 2 • 10 3 cm~ 2 sec fi. The rotation rates 
to achieve the line densities of order L = 10 6 • • • 10 8 
cm~ 2 , discussed below, are thus in the range 10 3 • • • 10 s 
rad/sec, i.e., of an order of magnitude reached in cen- 
trifuges only, so that this appears difficult to realize in 
experimental low temperature practice. In general, the 
turbulent vortex tangle will not be in a (local) thermo- 
dynamic equilibrium state. We may, however, expect 
that the contribution of line oscillations in the tangle to 
the internal energy will still be proportional to the line 
density Lq. Furthermore, any deviation from the pre- 
dicted dependence of the line oscillation contribution on 
temperature (for line densities which are not too high, 
so that the Kelvin spectrum remains valid), will give a 
measure of how far the turbulent tangle is away from 
thermodynamic equilibrium. 

The one-dimensional density of states for the Kelvin 
modes, within logarithmic accuracy, may be written in 
the form 
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For rings large compared to £ c , in the continuum ap- 
proximation of closely spaced excitation frequencies, this 
expression for the density of states then enables the calcu- 
lation of the specific heat of the Kelvin modes according 
to 
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There results the following expression for the dimension- 
less function in (Ef) 
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More explicitly, we may write for the specific heat of the 
Kelvin modes 
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the second line giving the absolute magnitude of the spe- 
cific heat expressed in terms of the temperatures, line 
densities, and sample dimensions of a typical experiment 
in *Re. There is no dependence on £ c (save for the very 
weak one contained in ln[n*]), as it should be for the large 
wavelength oscillation modes of a vortex line. 
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FIG. 1. Low temperature asymptotics of the total specific 
heat per volume V = V[^m 3 ], divided by VT = \/T[mK], 
for three different values of Lo/^/\n[n*] (the speed of sound 
c s = 240 m/sec). 



The ratio of the specific heats (|2|) and ([[]) is now given 
by 
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We may express this ratio in units of relevance for tur- 
bulence experiments in superfluid ^Re: 
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The contribution of line oscillations to the specific heat 
thus becomes noticeable, using the indicated line density 
of Lq = 10 6 cm -2 , for temperatures which are of the or- 
der mK. The total specific heat C = Cph + Ct per volume 
V, divided by y/T, is displayed in Fig. [1} The extrapola- 
tion of C I y/T to T = yields the smoothed line density, 
up to the very weak (square root of) logarithmic cutoff 
dependence on the largest possible number of waves on 
the filaments (the ultraviolet cutoff mode number), i.e., 
it gives the quantity Lq / y / Tnjn*J. 

In higher temperature superfluid turbulence experi- 
ments, values of L ~ 10 8 cm~ 2 have been reported |1. 
The number for Lq put into the estimates ([7J) and (g) 
above is thus rather conservative. The ultimate limit of 
achievable line densities in superfluid turbulence, and in 
particular the absolute line densities which are reached 
in turbulence experiments at mK temperatures, are as 
yet not clear H . Because the conventional second sound 
technique to detect vortices 0, fails at these tempera- 
tures, one has to look for different means to measure the 
vortex density. One conceivable possibility, using a plot 
like that shown in Fig. [l], is to measure the specific heat 
of the developed turbulent state in the mK range, di- 
vide by VT, and extract the (smoothed) vortex density 
from an extrapolation of the resulting asymptote to ab- 
solute zero. It should of course be mentioned that an 
experimental determination of specific heats as small as 
those estimated in equation (Q) is a less than trivial affair 
(though not entirely unrealistic). However, even if a di- 
rect determination of the specific heat of line oscillations 
Cl will prove to be difficult, this quantity is of relevance 
for the (thermo-) dynamical behaviour of the superfluid 
at the lowest temperatures. 

In the present context, it is also of interest to note that 
in the numerical simulation work of Nore et al. [^) , quite 
large line densities of order Lq ~ 10 10 • • • 10 11 cm -2 have 
been assumed. The resulting crossover temperatures, for 
which Cph = Cl, are consequently in these simulations 
already of order 100 mK. However, as already alluded to 
above, it is likely that for very large densities, modifica- 
tions of the Kelvin spectrum, by mutual induction of the 
vortex filaments, have to be taken into account. 

According to (||), the ratio of the contributions of 
Kelvin line oscillations and phonons to the specific heat 
scales with LqT~ 5 ^ 2 , so that simultaneously relatively 
low temperatures and high line densities are necessary to 
observe the influence of the specific heat of the Kelvin 
modes upon the thermodynamical behaviour of the tur- 
bulent superfluid. From equation (^J), we can conclude 
that the required values are not unreasonably large, and 



that the contribution of the Kelvin oscillations in the 
vortex tangle should indeed be observable. 
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